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Abstract 

A bound on the angle 7 of the unitarity triangle is derived using experimental 
information on the CP-averaged branching ratios for the rare decays —> tt^K^ 

and The theoretical description is cleaner than the Fleischer-Mannel 

analysis of the decays B^ —> tt^K^ and B^ —> vr^in that the two decay rates 
differ only in a single isospin amplitude, which has a simple structure in the SU(3) 
limit. As a consequence, electroweak penguin contributions and strong rescattering 
effects can be taken into account in a model-independent way. The resulting bound 
excludes values of cos 7 around ~ 0.6 and is thus largely complementary to indirect 
constraints derived from a global analysis of the unitarity triangle. 
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The study of CP violation in the weak decays of B mesons is the main target of 
present and future B factories. It will provide stringent tests of the flavor sector of the 
Standard Model and of the CKM paradigm, according to which all CP violation results 
from the presence of a single complex phase in the quark mixing matrix. The precise 
determination of the sides and angles of the unitarity triangle, which is a graphical 
representation of the unitarity relation V*i,Vud + Vcb^cd + ytl^td = 0 , plays a central role 
in this program [|l|. The angle [3 = —Pngiy^^Vtd) will be accessible at the hrst-generation 
B factories through the measurement of CP violation in the decay B —> J/'ipKs- The 
angle 7 = Arg(14*bf4d) is, however, much harder to determine. The sum {(3 + 7 ) can, in 
principle, be extracted in a theoretically clean way from measurements of CP violation 
in the decays i? —>• tttt (or in the related decays B ^ Tip and pp), but because of 
experimental difficulties such as the detection of the mode B —>■ this will be a 

long-term objective. 

Fleischer and Mannel have argued that some information on the angle 7 can be 
derived from the measurement of the branching ratios for the decays ti^K^ 

and B^ ti^K^, averaged over CP-conjugate modes [Q. In their original work, they 
obtained the bound 

r(5+) Br(5° ^ 7r-K+) + Bt{B^ tt+R-) .3 
t{B^) Br(i?+ —> 7 r+iF°) -|- Br(i?“ —>• 'k~K^) ^ ’ 


which would exclude a region around 7 = 90° provided that i? < 1, a possibility al¬ 
lowed by the hrst measurement of this ratio yielding R = 0.65 ± 0.40 |Q]. It has later 
been realized, however, that the bound o may be subject to sizable corrections arising 
from hnal-state interactions and electroweak penguin contributions, which are difficult 
to quantify in a model-independent way 0]-@]. To control such effects would require 
a more sophisticated analysis using information from other decays such as i? ^ KK 
1^, 1^. In addition to these theoretical obstacles, the prospects for deriving useful infor¬ 


mation on 7 using the ratio R are reduced by the fact that the CLEO Collaboration has 
announced a new measurement of this quantity yielding R = 1.0 ± 0.4 0.0 

In the present note we propose a variant of the Fleischer-Mannel analysis, which is 
theoretically cleaner and implies a non-trivial constraint on 7 provided present, prelim¬ 
inary data are con&med by future measurements. Electroweak penguin contributions 
play an important role in this analysis, but they can be controlled in a model-independent 
way. Consider the ratio 


Br(a-^ ^ + Br(a- ^ ^ . ,2 ,,, 

• 2|Br(B+ ^ n«K+) + Br(B- ^ ^ ' 

which would approach unity in the limit where b sqq penguin transitions involving 
top- or charm-quark loops dominate over the Cabibbo-suppressed h uus transitions. 
Deviations from this limit are measured by the quantity A*. The CLEO Collaboration 

^Some information on 7 can be obtained even when i? « 1 if one uses separate measurements of the 
decay rates for tt~K'^ and & K~ Q. 
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has recently reported the preliminary results —>■ = (1.4 ± 0.5 ± 0.2) x 10“^ 

and Bt(B^ = (1.5 ± 0.4 ± 0.3) x 10“® for the CP-averaged branching ratios 

|| 11 || , which imply 

= 0.47 ±0.24, AfP = 0.32 ±0.17, (3) 


corresponding to a deviation of i?* from unity of about two standard deviations. The 
goal of this note is to transform this measurement into a model-independent constraint 
on the angle 7 . 

The effective weak Hamiltonian governing the decays B —>■ ttK has the form [l^ 


n = 


^ f 10 

= ± E C.lKQ'i + KQi) - + h.c., 

V 2 t j=i^2 i=3 ^ 


( 4 ) 


where Xq = V^fjVqs are products of CKM matrix elements satisfying the unitarity relation 
An ± Ac ± At = 0, Ci are Wilson coefficients, and Qi are local four-quark operators. 
Relevant to the further discussion are the isospin quantum numbers of these operators. 
The current-current operators Q 12 ~ hsuu have components with AJ = 0 and A/ = 1; 
the current-current operators Qi 2 ~ hscc and the QCD penguin operators ~ 

have AJ = 0; the electroweak penguin operators Q 7 ,...,io ~ where Cq 

are the electric charges of the quarks, have AJ = 0 and AJ = 1. Since the initial B 
meson has 7=4 and the hnal states (vriP) can be decomposed into components with 
1=4 and 7 = |, the physical B —> nK decay amplitudes can be described in terms of 
three isospin amplitudes. They are called 7 ?i/ 2 , A 1 / 2 , and ^ 3/2 referring, respectively, to 
A7 = 0 with Bk = A7 = 1 with Bk = and A7 = 1 with Bk = i 


2 ’ 


2 ’ 


0 , 0- 


From the isospin decomposition of the effective Hamiltonian it is obvious which operator 
matrix elements and weak phases enter the various isospin amplitudes. The resulting 
expressions for the decay amplitudes relevant to our discussion are 


A{B+ ^ 7r+77°) = Ri/2 ± Ai/2 ± A 3/2 , 

-V2A{B+ ^ 7r°77+) = Bi/2 ± AI 1/2 - 2 H 3/2 , 

-A{B^ ^ 7r-77+) = B„2 - A„2 - A 3/2 . (5) 


Experimental data as well as theoretical expectations indicate that the amplitude 7 ?i/ 2 , 
which includes the contributions of the QCD penguin operators, is an order of magnitude 
larger than the amplitudes H 1/2 and ^ 3/2 P, 0- Yet, the fact that A 1/2 and ^ 3/2 are 
different from zero is responsible for the deviations of the ratios R and from unity. An 
important advantage of our analysis is that the two decay amplitudes entering the ratio 
R^ in (^ differ only in the sign of the isospin amplitude A 3 / 2 , which in the SU(3) limit 
is related to a matrix element of a single combination of local four-quark operators in 
the effective weak Hamiltonian. Despite the presence of contributions carrying different 
weak phases, there is thus a single strong-interaction phase associated with this ampli¬ 
tude. This fact will allow us to control rescattering effects and calculate electroweak 
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penguin contributions in a model-independent way. On the contrary, the decay ampli¬ 
tudes entering the ratio R considered in the original Fleischer-Mannel analysis differ in 
the signs of the two isospin amplitudes y4i/2 and ^ 43 / 2 , introducing hadronic uncertainties 
into the calculation @]-[§]. 

Taking advantage of the fact that the top- and charm-quark penguin contributions 
to i?i /2 are much larger than all other contributions to the isospin amplitudes, we write 
Bi /2 = \Pct\R'^^[R'^ RO{ey\, where is a strong-interaction phase, and is the weak 
phase associated with the top-quark penguin contribution. We generically denote by 
0{e) contributions arising from the quark decay b uus, which are Cabibbo-suppressed 
relative to the leading penguin amplitudes and carry the weak phase e*'^ of the parameter 
A„. Our main result is that the isospin amplitude A 3 / 2 , which causes the deviation of 
the ratio i?* from unity, can be written as 


3 ^ 3/2 


—^3/2 (e*'’' — (5 ew) , 


( 6 ) 


where is a strong-interaction phase, and ^ew is to a good approximation a real 

parameter accounting for the contributions of electroweak penguin operators. The fact 
that this parameter does not carry a non-trivial strong-interaction phase is crucial and 
will be derived below. We define the parameter £ 3/2 to be real and positive. In the 
diagrammatic amplitude approach of flavor-flow topologies, £ 3/2 = \T + C\/\Pct\ is the 
ratio of color-allowed plus color-suppressed tree amplitudes to the sum of the top- and 
charm-quark penguin amplitudes Using the parametrization (^, we obtain for 

the quantity A* defined in the result 


A* = -£3/2 cos A 0 (cosy - 5 ew) + 0{e '^), 


( 7 ) 


where A0 = 03/2 — 0p. The phase conventions adopted above take into account the 
naive phase relation between penguin and tree contributions, such that all non-trivial 
strong-interaction phases are accounted for by the phase difference A0. 

The value of the parameter £ 3/2 can be estimated using various experimental and 
theoretical information. The penguin contribution \Pct\ can be extracted from the CP- 
averaged branching ratio for the decays —>■ which are expected to receive 

only a very small contamination from other contributions such as up-quark penguins 
or annihilation diagrams. Using the experimental value for this branching ratio quoted 
earlier yields \Pct\ = (3.74±0.72) x 10“^. (We give amplitudes in “branching ratio units”, 
so that Br = |Alp.) The tree contribution |T -|- Cl can be estimated experimentally 
invoking SU(3) symmetry, in which case it can be extracted from the branching ratio for 
the decay 


7r^7r° 


using the relation 


T + C = -V 2 ^ ^ A{B+ 7r+7r°), 

^ud Jtt 


( 8 ) 


where the factor fx/fn = 1-22 ± 0.01 accounts for the leading (i.e., factorizable) SU(3)- 
breaking corrections. The CLEO Collaboration has recently reported a signal for this 
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decay mode with a significance of 2.3 standard deviations. They quote the upper limit 
Br(i?=*" — >• 7r^7r°) < 1.6 x 10“® at 90% confidence level [|TI|, which implies \T + C\ < 
1.56 X 10“^. Taking this signal seriously, we may derive from the reported event rate and 
efficiency the branching ratio Br(i?^ —>■ = (0.59lo!27) ^ 10 “^, which yields the value 

\T + C\ = (0.94]'^Q;2f) X 10“^. The tree contribution may also be obtained theoretically 
employing the factorization hypothesis. Using the values of the hadronic form factors 
determined in Ref. combined with the new world average \Vub\ = (3.56±0.56) x 10“^ 


we get |T + U| = (0.88±0.23) x 10 in good agreement with the estimate obtained 


using experimental data. We take the average of the two numbers and combine it with the 
experimentally determined value for the penguin amplitude to obtain £ 3/2 = 0.24 ±0.06. 
It is important that, once the branching ratio for the decays —> vr^vr^ and 

are measured more precisely, the parameter £ 3/2 can be determined with smaller 
errors and without almost any recourse to theory. Note that with the above value of 
£ 3/2 the terms of 0{e^) omitted in ([^ are expected to be reasonably small compared 
with the sizable experimental uncertainties in (^). Nevertheless, we shall comment on 
the relevance of such higher-order terms below. 

Inserting the value of £ 3/2 into relation (|^) and neglecting higher-order terms, we 
obtain 


cos A 0 (5ew ~ cosy) = 


^3/2 


= 1.33 ±0.78 


|<5ew ~ cosyl > 1.33 ± 0.78 . 


( 9 ) 


The last step missing to turn this result into a useful constraint on y is an estimate of 
the contribution from electroweak penguins to the isospin amplitude % 3 / 2 - This can be 
done without encountering hadronic uncertainties. The crucial observation is that the 
electroweak penguin operators Qg and Qio, whose Wilson coefficients are enhanced by 
the large mass of the top quark, are Fierz-equivalent to the current-current operators 
Qi and Q 2 - As a result, the AJ = 1 part of the effective weak Hamiltonian for B —> iiK 
decays can be written as p, 

nAi=i = ^ { (a.C'i - ^AiCg) Qi ± (a„C '2 - ^AiUio) Q 2 ± .. j ± h.c., (10) 

where Qi = \{Q^ — Qf). The dots represent the contributions from the electroweak 
penguin operators Q 7 and Qs, which have a different Dirac structure. The Wilson 
coefficients of these operators are so small that their contributions can be safely neglected. 
It is important in this context that for heavy mesons the matrix elements of four-quark 
operators with Dirac structure (U —A) 0 iV ± A) are not enhanced with respect to those 
of operators with the usual (U —A)®(U—A) structure, as can be checked explicitly in the 
factorization approximation. To an excellent approximation, the net effect of electroweak 
penguin contributions to the AJ = 1 isospin amplitudes m B ^ ttK decays thus consists 
of the replacements of the Wilson coefficients Ci and C 2 of the current-current operators 
with the combinations shown in (^). In the SU(3) limit, [/-spin invariance implies that 
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the isospin amplitude A 3/2 is given by the matrix element of the combination (Qi + Q2) 
only, the coefficient of which is |[A„(Ci+ C 2 ) — |Ai(C '9 + C'io)]. The difference, {Q1—Q2), 
does not contribute to this amplitude. To see this, consider the relation (||) between the 
current-current part of the amplitude ^ 3/2 and the amplitude for the decay 
Bose symmetry requires the two pions to be in an / = 2 state, so that only the AI = | 
part of the effective weak Hamiltonian for B —> tttt decays contributes to the decay. Since 
the combination {Q[ — Q 2 ), where the prime on the operators indicates the replacement 
s —>■ d appropriate for AS = 0 transitions, has AI = A, the decay B~^ —>■ 7r~^7r^ receives 
a contribution from the combination (Q[ + Q 2 ) only. SU(3) symmetry then implies that 
it is only the combination (Qi + Q2) which contributes to the isospin amplitude ^ 3/2 in 
B — tiK decays. It follows that the quantity ^ew in (H) is given by 


^EW — 


3 Cg + Cig 

2xmb Cl + ^2 


0.66 ± 0 . 11 , 


( 11 ) 


where we have used A„/Ai ~ with A ~ 0.22 and Rf, = X~^\Vub/Vcb\ ~ 0.41 ± 

0.07. The ratio of the Wilson coefficients is to a very good approximation independent 
of the choice of the renormalization scale. We take {Cg + Cio)/(C'i + C 2 ) = —1.14a, 
which is obtained using the leading-order coefficients^] at the scale fi = rrib [j^j. For 
the electromagnetic coupling at the weak scale we take a = 1/129. Note that the main 
uncertainty in the value of (5 ew results from the uncertainty in Vub, which is likely to be 
reduced in the near future. 

SU(3)-breaking corrections to the relation O are controlled by the ratio of the 
following operator matrix elements: 


{irA'(/ = i)IQi+Q2|B+) 


( 12 ) 


The magnitude of this ratio can be estimated using the generalized factorization hypoth- 
p0[| to calculate the matrix elements of the current-current operators. This leads 


esis 

to 


<5su(s) = I -7 « 1 - 3 %, Av = 0, (13) 

where ax = fK{fn% — {m\) and = / 7 r(m| — {rri^) are com¬ 

binations of hadronic matrix elements, and ~ 0.45 is a parameter controlling non- 
factorizable corrections. Despite the fact that the factorization approximation may be 
crude, this estimate suggests that SU(3)-breaking corrections to the result m are very 
small. To linear order in SU(3) breaking, the results (j^) and (® can be corrected by 
replacing ^ew with the effective value 


iceff f, u, cos{A(j) + Aip)\ 

^EW - I 1 - fc asu( 3 )- cosA(j) - ) 

^Using next-to-leading order coefficients would change this value by less than 2%; however, at this 
order the coefficient functions may depend on the choice of the operator basis, so one has to be careful 
when using Fierz indentities. 
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Table 1: Constraints on 7 for different valnes of the ratio A*/£ 3/2 


A*/£ 3/2 

deviation from central value 

bound on cosy 

excluded region 

2.11 

Tier 

no solution 

all 

1.33 

central value 

< -0.55 

|y| < 123° 

0.55 

—la 

< 0.23 

|y| < 77° 

0.16 

—1.5cr 

> 0.64 or < 0.62 

50° < |y| < 52° 


where 


, Ca - Cl Cg - Cio 
k = ^ ^ ^ 3.43 


(15) 


Ci + Ci C9 + C10 

is a combination of Wilson coefficients. Since {Q 1 — Q 2 ) and {Q 1 + Q 2 ) are local operators, 
it is unlikely that the strong-interaction phase difference A(p could acquire any significant 
value. We thus expect that the main effect of SU(3) breaking is to reduce the value of 
( 5 ew in (|Tll) by an amount of order 5%. To account for this effect, we shall from now 
on use the value ^ew = 0.63 ± 0.15 with an increased error, which is large enough to 
cover possible small contributions from a non-zero phase difference A(p or deviations 
from the factorization approximation. It is interesting to note that our general results 
for the structure of the electroweak penguin contributions to the isospin amplitude A 3 / 2 , 
including the pattern of SU(3)-breaking effects, are in full accord with a model estimate 
by Deshpande and He However, these authors did not realize that such contributions 
can be calculated in an essentially model-independent way. As a result, electroweak 
penguin effects can be taken into account in the method proposed in Ref. for learning 
7 by comparing the rates for the processes —>■ ^ and 7 r''' 7 r‘^ with their 

charge conjugates. This will be discussed in more detail elsewhere . 

The role of electroweak penguins in the bound @ is an important one, and it is 
crucial that the parameter ^ew can be calculated in a reliable way in terms of perturbative 
Wilson coefficients and measured quantities. The result is that the constraint (j^) excludes 
values of I 7 I around arccos(0.63) ~ 51°. More precisely, the excluded region depends 
on the value of the experimental number for the ratio A*/£ 3/2 on the right-hand side 
of the bound. The allowed regions for cosy must satisfy cosy > (5 ew + ^*/£z /2 or 
cosy < (5 ew — A*/£ 3 / 2 . In evaluating these results we lower or increase the value of 
<5 ew by one standard deviation so as to be conservative, i.e., we use (5 ew = 0.48 in 
the hrst relation and (5 ew = 0.78 in the second one. The hrst solution only exists if 
A*/£3/2 < 1 — 5 ew- If A*/£3/2 > 1 + ( 5 ew there is no solution for y at all. The results 
obtained for some representative values of A*/e 3/2 are shown in Table 2. It is remarkable 
that the constraints on y derived from the bound (P) are to a large extent complementary 
to the preferred region for this angle obtained from a global analysis of the unitarity 
triangle, using information from semileptonic B decays, B-B mixing, and CP violation in 
the kaon system. A typical range of values allowed by such an analysis is 47° < y < 105° 
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PI 


where the upper bound is determined by the improved lower limit on the quantity 

For the central value of A*/e 3/2 our bound is 


Arris controlling Bg-Bg mixing [H 


inconsistent with this region, while for the la lower value 82,12 = 0.55 it would 
exclude values of I 7 I below 77°, thus leaving the rather narrow range 77° < 7 < 105°. 
Clearly, a more precise measurement of the parameter A*/£ 3/2 could potentially yield a 
very non-trivial constraint on 7 and provide a stringent test of the CKM paradigm. 

The weakest point of our analysis appears to be the neglect of higher-order terms 
in £ in the expression for the quantity A* in (|^), given that £ 3/2 ~ 0.24 is not such a 
small parameter. We recall that £ is a generic measure of the strength of the Cabibbo- 
suppressed b —>■ uus transitions, which carry the weak phase e*'’', relative to the leading 
penguin transitions. Whereas a full account of all higher-order terms of this kind would 
introduce unknown strong-interaction parameters, we now show that the potentially most 
important terms can be analysed without much additional complication. The reason 
is that contributions proportional to the weak phase in the decay amplitude for 
B^ —> are likely to be much smaller than those in the amplitude for B'^ —> , 

because they could only be induced via hnal-state rescattering through soft annihilation 
or up-quark penguin topologies 0. Model estimates indicate that those effects typically 
lead to contributions of order a few percent [0-|@]. In other words, if we write i?i /2 + 
A 1/2 + A 3/2 = we expect that Sa -C £ 3 / 2 - This observation justifies 

treating these rescattering contributions from a numerical point of view as quadratic 
rather than linear in £ 3 / 2 . The generalization of the relation (|^) to the next order in £ 3/2 
then becomes 


A* = £3/2 cos A 0 (^Ew - cos 7 ) 


'3/2 


(3cos^A0 - 1)(5ew - C 0 S 7 )^ - sin% -h 0(£3/2, £a£ 3 / 2 ) • (16) 


The question is whether the terms of order £ 3^2 could increase the value of A*, thereby 
weakening the bound on 7 derived in this letter. Since £3/2 -C 1, the quantity A* takes 
its maximum value for | cos A^l = 1 , so that 


A* < £ 3/2 |5ew - cos 7 1 



(17) 


If experimentally it is found that A* > £ 3 / 2 / 7 ^ 0.17, the second term is negative and 

strengthens the bound on 7 . Provided future measurements confirm the present value of 
A* in (|^) within one standard deviation, the bound (|^) derived in linear approximation 
is thus a conservative result. On the other hand, if A* would turn out to be much smaller 
than the current central value, the bound may be weakened by higher-order terms, but 
those would generally be smaller than £ 3 / 2/2 3%, which is an almost negligible amount. 

In summary, we have shown that a non-trivial bound on the angle 7 of the unitar- 
ity triangle can be derived from a measurement of the branching ratios for the decays 
B^ —>• ttK, averaged over CP-conjugate modes. Electroweak penguin contributions. 
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which play an important role in these decays, can be controlled in a model-independent 
way using Fierz identities and SU(3) symmetry relations, with strong indications that 
SU(3)-breaking corrections are very small. 
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